Anomalous thermal expansion of SboTe.^ topological insulator 
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We have investigated the temperature dependence of the linear thermal expansion along the hexagonal c axis 
(AL), in-plane resistivity (p), and specific heat (C p ) of the topological insulator Sb 2 Te3 single crystal. AL 
exhibits a clear anomaly in the temperature region 204-236 K. The coefficient of linear thermal expansion a 
decreases rapidly above 204 K, passes through a deep minimum at around 225 K and then increases abruptly 
in the region 225-236 K. a is negative in the interval 221-228 K. The temperature dependence of both a and 
C p can be described well by the Debye model from 2 to 290 K, excluding the region around the anomaly in 
a. 



Topological Insulators (TIs), being a topic with 
fascinating physics and high application potential, have 
been a highly emerging field in last five yearaM These 
materials possess an insulating energy gap in the bulk 
due to the spin-orbit coupling and gapless surface states 
protected by the time-reversal symmetry. The topologi- 
cal^ protected surface states host a variety of properties 
like robust conductance, relativistic Dirac cone, spin- 
momentum chirality, and quantized anomalous Hall 
effectPEl. These dissipationless surface states are poten- 
tially useful for applications in spintronics and quantum 
computations, and also in exploring the fundamental 
physics like Majorana Fermions, magnetic monopoles 
and axion electrodynamics. Though, many members 
of TIs such as Bi 2 Se 3 , Bi 2 Te 3 , Sb 2 Te3, etc. have been 
studied earlier as good thermoelectric materials, it is 
essential to characterize their basic physical properties 
more accurately down to the lowest temperature for 
device designing. The atomic arrangement of TIs is 
a layered structure, with each layer consisting of five 
monoatomic planes of Te(Se)-Sb(Bi)-Te(Se)-Sb(Bi)- 
Te(Se). Each of these five atomic planes is referred to as 
a quintuple layer (QL). The lattice properties of these 
materials are of considerable interest because of the 
weak van der Waals interaction between the adjacent 
QLs. Previous reports show that the thermal expansion 
in these materials is quite anisotropic and exhibits 
an anomalous temperature dependence, presumably 
because of the weak bonding 3 6 . Chen et studied 
the temperature dependence of the lattice parameters 
of Bi 2 Se3 and Sb 2 Te3 single crystals. Below 150 K, the 
extracted coefficient of linear thermal expansion a can 
be described well by the Debye model. While above 150 
K, the estimated value of a is larger than the theoretical 
one and does not show saturation-like behavior well 
above the Debye temperature. Furthermore, from the 
Raman Spectroscopy measurements, it has been shown 
that the thermal expansion contribution term Acj(T) 
accounted for 40% of the total phonon frequency change 
with temperature^. 

In most cases, as the a was determined from the 
temperature variation of the lattice parameters, the 



magnitude and nature of the T dependence of the re- 
ported anomalies differ significantly from each other due 
to the lack of accuracy. Thus an accurate measurement 
of thermal expansion is necessary to shed some more 
light on the nature of anomaly and see whether it has 
any correlation with electronic properties. In this regard, 
we have measured the temperature dependence of the 
in-plane resistivity (p), specific heat (C p ) and linear 
thermal expansion (AL) along c axis of Sb 2 Te3 single 
crystal. We observe that AL(T) exhibits an anomaly in 
the region 204-236 K. No anomaly is observed in C p (T) 
and p(T) data. Both C„(T) and a(T) can be described 
well by the Debye model. 

Single crystals of Sb 2 Te3 were grown by melting the 
stoichiometric mixture of elemental Sb (99.999%) and Te 
(99.999%) at 850 °C for 24 h in a sealed vacuum quartz 
tube. The sample was cooled over a period of 48 h to 620 
°C and kept at this temperature for 96 h before quench- 
ing in liquid nitrogen. Crystals of millimeter size with 
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FIG. 1. Temperature dependence of the in-plane resistivity 
p of the St>2Te3 single crystal. Inset shows the power-law fit, 
p(T) — po + aT n with n~4/3, to the resistivity data in the 
temperature region 20-270 K. 
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shiny flat surface can easily be obtained from the boule 
of Sb2Te3. The phase purity of the sample was checked 
by powder x-ray diffraction (XRD) method with CuK Q 
radiation in a Rigaku x-ray diffractometer (TTRAX If) 
and no impurity phases were detected. The XRD pattern 
can be fitted well with a hexagonal structure with space 
group Rim. The evaluated lattice parameters are ah ex 
= 4.2654 A and c^ ex = 30.4438 A. Four probe measure- 
ments of the electrical resistivity between 2 and 300 K 
were done using a commercial cryostat (Cryogenic Ltd.). 
For the thermal expansion measurement, a miniature 
tilted plate capacitance dilatometer was utilized which 
allows an accurate study of sample length change^. We 
have measured the macroscopic length change AL(T) of 
a crystal of ~1 mm thickness and calculated the coeffi- 
cient of linear thermal expansion, 



J_ d{AL) 

L ' dr ' 



(1) 



where Lq is the length of the sample at lowest measured 
temperature. The specific heat was measured by the re- 
laxation time method using a Quantum Design physical 
property measurement system. 

Figure 1 shows the temperature dependence of p for 
Sb2Te3 single crystal. With decreasing T, p decreases 
continuously down to 20 K and then saturates to a con- 
stant value. No anomaly has been observed in the mea- 
sured temperature range. Both the magnitude and T 
dependence of p for the present sample are comparable 
to that for Bi2Se3 and Bi2Te3 crystals with carrier den- 



sity 10 i8 -I0 ly cm 
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In Bi2Sc3, for carrier density 



greater than 10 18 cm" 3 , p(T) data reflect good metallic 
behavior^. In comparison to the above mentioned crys- 
tals, the residual resistivity ratio /?3ook/P2K ^5 for the 
Sb2Te3 crystal is higher 9 " 11 . In the temperature region 
20-270 K, the resistivity can be fitted well by a power- 
law expression, p(T)=po + oT n with exponent n~ 4/3 as 
shown in the inset of Fig. 1. The observed value of the 
exponent does not seem to classify this compound into 
hither-to-known metallic systems. Normally, in metal- 
lic elements, alloys, and compounds, p shows T 2 depen- 
dence at low temperature and T-linear behavior at high 
temperature. In magnetic materials, p may show T 3 / 2 
dependence at low temperature due to the magnetic con- 
tribution to the scattering. In the bulk samples of TIs, 
the transport properties of the surface states are often 
mixed with those from the bulk states making it difficult 
to observe insulating behavior in p(T). As the bulk con- 
ductivity decreases whereas the conductivity due to the 
surface states increases with decreasing temperature, the 
electrical conduction at low temperature is dominated by 
the surface. 

Figure 2(a) shows the temperature dependence of 
AL/Lq along the hexagonal c axis of Sb2Te3 crystal in 
the range 4 to 290 K. We have also measured the ther- 
mal expansion in presence of magnetic field but AL/Lq 
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FIG. 2. (a) Temperature dependence of the linear thermal 
expansion AL/Lo along the c axis of the Sb2Te3 single crystal. 
Inset shows the expanded view of the anomalous region, (b) 
Temperature dependence of the coefficient of linear thermal 
expansion a. The solid line is the fit to the Debye model 
[Eqn.(2)]. 



does not show any H dependence up to 7 T. The length 
decreases with decrease in T as it occurs in most met- 
als, however, AL/Lq exhibits a broad minimum at ^228 
K followed by a broad maximum at ~221 K [inset of 
Fig. 2(a)]. With further decrease in T, AL/Lq decreases 
continuously down to the lowest measured temperature. 
This indicates that there exists a small temperature re- 
gion 221-228 K where an anomalous thermal expansion 
occurs. The anomalous behavior is clearly reflected in 
the coefficient of linear thermal expansion. In Fig. 2(b), 
a is shown as a function of T. Except for 204 K<T<236 
K, a increases smoothly with the increase of T and sat- 
urates at high temperature. Above 204 K, a decreases 
sharply with increasing T, passes through a deep mini- 
mum at around 225 K and then increases abruptly in the 
region 225-236 K. For 221 K<T<228 K, a is negative. 
Though the anomaly in the region 204-236 K is slightly 
asymmetric, it is not a A-like as in the case of a second- 
order phase transition. Excluding the region around the 
anomaly, the a(T) curve can be fitted well in the whole 
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FIG. 3. Temperature dependence of the specific heat of the 
Sb2Te3 single crystal. The solid line is the fit to the Debye 
model (see the text). Inset: C p /T data below 5 K has been 
fitted (solid line) with the expression C p /T = 7 + A3T 2 + 
A 5 T 4 . 



over a wide range of T above 350 K and observed a 
very weak T dependence in C p . The specific heat data 
in Fig. 3 can be described well by the Debye model 
(Eqn. 2) where evo is replaced by 9NR and an extra 
term jT due to electronic contribution is added, where 
7 is the Sommerfeld coefficient. Also, the estimated 
value of 0£> ~165 K agrees very well with that reported 
in literature^. For the accurate determination of 7, 
the specific heat data below 5 K can be described by, 
C p = C ei + C ph = 7T + A 3 T 3 + A 5 T 5 , where A 3 and A 5 
are the coefficients of the phononic contribution (inset 
of Fig. 3)P2I. From the fit, the value of 7 can be confined 
to 0.7±0.7 mJ mob 1 K~ 2 . 

In conclusion, we have investigated the temperature 
dependence of the resistivity, specific heat, and linear 
thermal expansion of the topological insulator Sb2Te3. 
AL exhibits an anomalous thermal expansion in the 
temperature region 221-228 K. p(T) shows an unusual 
power-law behavior for 20 K<T<270 K. Both a(T) 
and C P (T) can be described well by the Debye model 
excluding the region where anomaly is observed in a. 



temperature region with the Debye model, 

where ag is T-independent fitting parameter and 0£> is 
the Debye temperature. The estimated value of Od is 
235 K. Chen et alW determined a from the temperature 
dependence of lattice parameters of Sb2Te3. They ob- 
served that a(T) can be described well by the Debye 
model up to 150 K with 0£>~2OO K. However, above 150 
K, a along c axis initially decreases slowly with T down 
to 180 K and then increases, as a result, the experimen- 
tal value of a is higher than the Debye model for T > 
230 K. On the contrary, a perpendicular to the c axis 
increases monotonically with T. This suggests that the 
large anomaly in thermal expansion along the c axis is 
associated with the inter-layer van der Waals interaction. 

In order to investigate whether this anomaly in 
thermal expansion is due to the phase transition, we 
have measured the specific heat of Sb2Te3 from 2 to 300 
K as shown in Fig. 3. The C P (T) data do not show 
any anomaly around the temperature interval where 
the anomalous thermal expansion is observed, ruling 
out the possibility of any electronic or structural phase 
transition. Also, one can see that C p tends to saturate at 
~125 J mol -1 K _1 at high temperature which is close to 
the Dulong-Petit value. Both the value and temperature 
dependence of C' p for the present sample match very 
well with the reported specific heat data for Sb2Tejpl 
Pashinkin et aZ™ measured C p on Sb2Te3 single crystal 
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